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1. Int reduction . In a Ikniard-t ype convert Ion problem mu* seeks, 
e,g», t» determine tin* nt: at. tonary flows of an infinite Inver of lluid 
Ivin? 1 , between two rigid horizontal walls and heated uniformly from below. 
Sin*!i a problem possesses a unique, motionless conduction solution when 
the parameters of the problem lie within a certain range but, as the 
temperature difference across the layer increases beyond a certain value, 
other, convective, motions appear. These motions are often cellular in 
diameter In that their streamlines are confined to certain well-defined 
"cells" having, e.g., the shape of rolls or hexagons. The purpose of 
(his paper is to formulate a "selection principle" that explains why 
hexagonal cells seem to be "preferred" for certain ranges of the parameters, 

Uenard-type problems and their generalizations play an important role 
In fluid dynamics and have been investigated in recent years by a number 
of authors. Convection problems have been studied, e.g., by Sehluter, 
l.ortx and Busse f 19 1 and Fife and Joseph [A] using expansion methods, by 
Busse fl] using variational methods, by Kirehgassner [9,10] using the 
Ivapunov-Schniidt method, by Sattinger [ 17 , 18 ] and Golubitskv, Swift and 
Knob loch [5] using group-theoretic methods, and by Buzano and Golubitsky 
| .’ j using group-theoretic methods and singularity theory. The reader 
Ik referred to the above papers and to the book of Joseph [ 7 ] for a com- 
prehensive introduction to Benard-type problems. 

An important aspect of the work of Busse [1] is that the "extremum 
principle" and the stability results there are independent of the number 
of critical wave vectors corresponding to a given critical wave number. 

In the same spirit an important aspect of this work is the formulation 
and verification of a selection principle in a setting that is independent 
of any fixed number of critical wave vectors. Although our study is 
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rest rioted to functions doubly periodic In the horizontal plane, the 
(finite) number of critical wave vectors can be taken arbitrarily largo 
hv proper choice of the period rectangle. Moreover, In the case of tin.' 
hexagon; 1 lattice this choice can be made in .such a wav that the critical 
wave number and the "size" of the resulting hexagonal cells are kept 
lixed. Thus, whereas other methods offer a complete bifurcation analysis 
mi t lie hexagonal lattice in the usual six-dimensional setting, Lhe methods 
of this paper prove useful for a stability analysis on the hexagonal 
lattice in the general case of an arbitrarily large number of critical 
wave vectors (see also the discussion in Section 7) . 

To obtain a physical interpretation of the extremum principle i ■ 

Ml, Palm [15] derived in the time-dependent problem a minimum principle 
tor a type of generalized dissipation, V, namelv that, as time increases, 

V decreases and attains a minimum value on r >y state solutions (see 
I 15, p. 7.4 1 4 | ) . To treat tine generalized Uenard problem studied here, 
we introduce an analogous sort of functional, called the generalized 

dissipation (see (3.23) in Section 3 below). It can be shown for time- 
dependent problems in a formal way as in [15] that the associated time- 
dependent V decreases as time increases and assumes a minimum on steadv 
state solutions. Since 1/ = 0 for the motionless conduction solution 
and since \J initially increases in the steady state problem along a 
siiberitical branch of convective solutions bifurcating from the conduction 
solution at the critical Rayleigh number, R c> it is natural to conjecture 
L hat: what we shall call a "selection principle" is related to the existence 
of a convective solution for which 1/ - 0. Presumably, such a solution 


would correspond to a point on an "upper" branch because V > 0 on "lower" 
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subcritical branches. Using such an interpretation, one could replace 
the formal "geometrical" condition for upper branches used in [1, p . 6 3 1 ] 
b v the exact analytical condition V * 0. This would be an important 
l list .step in n inability analysis since subcrilical solutions Iving 
on upper branches are the ones most .likely to be stable. 

The basic idea of the paper can now be stated as follows (see also 
the related but somewhat easier approach used in [12] to solve a class 
oi variational problems arising in nonlinear shell theory — the parameter 
; in (,!.')) plays the role of the "structure" parameter « in [121). 
Instead of solving only the Boussinesq-type equations given in (2.1) as 
is usually done, we solve the equations in (2.1) together witli the eon- 
4 rain! that, 'or fixed t near y - 0, 1/ — 0 is a local minimum of 
l\ One anticipates here that the condition l'* = 0 will lead to a solu- 
tion on an upper branch and that the minimization condition will lead 
to a stable solution. In this paper we show that such an approach does, 
in tact, yield stable, subcritical solutions of the generalized Benard 
problem, when Y is sufficiently small. Such solutions may even be 
considered as "large" solutions because they are both subcritical and 
stable whereas "small" subcritical solutions bifurcating from the con- 
duction solution at R £i are always unstable. In this sense our method 
mnv he regarded as a "selection principle" for obtaining "large", stable, 
subcritical solutions because the. method selects certain solutions of 
equations (2.1) while excluding certain others. By "stability" here and 
throughout the remainder of the paper we mean "linearized stability" 
relative to some appropriate Hilbert space. 

The outline of the paper is as follows. In Section 2 we give an 
operator-theoretic formulation of a certain tvpo of generalized Benard 
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problem and in Sections i d 4 we reduce the given infinite- dimensional 
problem to one of Halving a finite-dimensional system of equations, the 
so-called He 1 eel ton equation:- The selection equations are derived 
bv means of splitting techniques such as those used in the Lyapunov- 
Schmidt meLhod in bifurcation theory but the equations obtained are not 
the usual bifurcation equations associated with the problem. The works 
of Kirchgassner [10] and Sattinger [17] play an important role in these 
preliminary sections. Sections 5 and 6 contain the main results of the 
paper. In Section 1 we solve the selection equations in a general 
setting by the use of variational methods and present a linearized stabilit 
analysis of the resultant stationary flows. In Section 6 we show for 
the hexagonal lattice that the classical hexagonal cellular solutions 
are generated from the absolute minimum of an appropriate selection 
functional and that such a minimization property is independent of the 
dimension of the basic underlying finite-dimensional problem. Thus, 
since the classical hexagonal cellular solutions are also stable, they 
are in some sense the preferred subcritical convection solutions. 


OF POOil 


v 


* J 

n* 

4 


> . 

2 , Formulation of the problem. In this section we formulate n 
general J />■ I Ildnard problem for certain temperature- dependent fluids .uni 
introduce a Hilbert, space setting for its study. Tin* particular problem 
described below Is chosen mainly lor convenience. The methods of the paper 
apply also to a much wider class of convection problems (e.g., see [1]), 

The generalized Bernard problem studied here is to determine the 
stationary flows or an infinite layer of fluid between two rigid, horizontal 
walls, and heated uniformly from below, The fluid density, r, is assumed 
to he constant, say > = p^, except in the gravity term where it is taken 
to he quadratic in the temperature, T, i.e., 


s) = e o rt - a (T - T q ) - b (T - T 0 ) 2 ], 

wln re T f) Is the average of the (constant) temperatures T^ on the upper wall 
ami Tj on the lower wall. Under this assumption on f, one is led, after 
scaling the variables suitably, to the system of Bousslnesq-type equations 
given in (2.0 below. The equations relate, at each point of the set 


W e fx » (x,y,z) < x,y ■ 



the fluid velocity vector, u = (u^,U 2 ,Ug), scalar pressure, p, 
variable, 0, measuring the change in temperature from its value 
conduction state (see, e.g., [9] where u, p, 0 are. related by 
lu those used here ): 


(2.1) (a) -Au - Akf^(O) + Vp = -(u*V)u + kf^O) 

(b) -(Pr) _1 A0 - Au 3 = -u • VO 

(c) V • u = 0 

(d) u = 0, 0 = 0 for z - ±-|- . 


and the sc lar 
for the pure 
a factor ' 
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In (2.1), k * (0,0,1) , V * f’.*^ , Vy , j, and A is tin* Lap 1 not* operator; tin* 

Pram! t 1 number, Pr, equals the ratio of kinematic* viscosity, *, to thermal 

minluai ivity and is regarded ns a fixed constant t hro ug! unit the paper; the nr is- 

o 

hoi i number Or * agd(T,j - TgJ/v" (g * gravitational constant, cl r - thickness 
ol i lie unsealed layer); A * v^cTr and 

(2,2) (a) f j (* ) = '* ( 1 - 2, s) 

(b) f ? (*0 * r' ;2 , 


where i is a "structure" parameter given by 
(2. U Y * b (T 1 - T 2 )/a, 


the Itav'li* i }>h number, Rn, Is related to A by Ra - PrOr « Pr> , 

We shall seek solutions having a doubly periodic cellular structure*. Thus, 
riven positive numbers and ^ (to be specified below), we set 


I* - (x * (x,y,z):0 < 


1 


0 < y < 


2 r 


1 


. 1 


1 


We next introduce the (complex) Hilbert space, H, defined as the closure of the 


set 1 v - (u^ ,u., ,u^,o) : v smooth, periodic in x with period “-* , periodic 

•j ^1 

in v with period “ , v - 0 in a neighborhood of \/.\ ~ , j and ' >v 0 1 in 

1 ’ •> 

the iiorii! 1| • |j associated with the inner product 


(v,w) = 


3 

[ T. Vv . • Vw . 

h j=r J " J 



Here and throughout the paper a bar over a quantity denotes complex conjugation 


am! tlie symbol 


111 1 . )ll 2 . S 

+ Ty- + ST 


: i. .1 JL n 

hx 5 <ly 5 Hz’ j 


when used with elements of H. Thus 


It we take the scalar product of (2.1n,b) with w < H, use (2.1e,d) and 
i n l egrat i on by parts, then for v = (u,h) we obtain 


( 2 . i ) (v,w) - A (L v,w) = (F. (v),w). 

y r 

Here the linear operator ly’-ll v H and quadratic operator F^tf/ -> H are given by 
L, y = L - r'M , 


( 2 . 2 ) 


F » v + >(’. 
Y 




(2.7) (Mv,w) * 2^-v^w^ 

t.’.K) (F(v) ,w) ~-j (v*Vv)*w 

C.'.'J) (H(v) ,w) ~ ( V/ ) 2 w.. 

for .ill v,w in H. Since in (2,4) w is an arbitrary element of H we ^ 
that a smooth solution v * (u,0) of (2.1) in H satisfies tilt* operator equatien 

( * ) 0 * v - XL.^v - F (v) , v H, X < 1R^ , 

In fact, one can apply standard regularity methods (e . g. , see [11,11,14]) Lo 
show that problems (2.1) and (*) are equivalent. 

In order tro study solutions of (*) wo shall require properties of the 
linearized version of (*) when 1=0, 

C\ 10) 0 « v - Xr.v, v < Ilf X . B 1 . 

Tlie linear eigenvalue problem (2,10) :,s equivalent to the classical problem, ter 

2 -[I 2 T 

smoolh u, p, 0 periodic with periods ~~ in x and - w in y, obtained bv 

1 '2 

omitting the nonlinear terms in (2.1), This linear problem is well studied (sut*,e.c., 

[5.o,I0,ll|). The eigenfunctions are complete in H and are obtained from the 
relations k » (k 1 ,k 2 ,0), a = (k 2 + k 2 ) 1/2 , i * A~T and 


(2.11) (a) Uj = e ik A? 4).(z), j = 1,2,3, 

\ ,, ik*x . 

(h) 0 = e •• rj« 4 (5s) , 

, \ ik*x -22,, 

(o) p ~ e - - a D tf»3 , 

(d) <!>. = ia” 2 k^^, 


(:) = 1 , 2 ) . 
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•) $ 9 (I 

I> 9 - 5% a prime denotes r nnd and * satisfy 

dz u/ } * 4 



IP) 

(a) 

0 B 

»% - 

j 

A * . 

' 4* 




(b) 

0 * 

1^4 

+ H.j, 




(e) 

*1 ■ 

■ !, 1 * ,1 
•3 *2 

K 

O 

S3 

rr 



Out* ran allow (o.g. , sot* [6]) for m > 0 that the eigenvalue problem 12,12) 

h.ei a countable .sequence of positive, simple eigenvalues, 0 ^ 

depend Ing continuously on M, Moreover, (>7) ■* '« as M M)* or - •» , r.«n- 

oequenl ly, U-j (*•) assumes an absolute minimum at some ^ > 0 depending only 

• 'll the Prandtl number, Pr. We as sume t hrough the paper that o ts unique 

1 b n n ( ( - ) ’• if * J 4 M () . (This property is suggested by numerie.il t 1 1 — 

cul itions | 5J and is usually assumed for Uenard-Lype problems.) lor clv.n in- 

I eg,e rs n^, nip we now choose ■>., such that 

r , , n .2 2 2 , 2, 2 
1 l) * n 0'| + "'O'' 0 ’ 

In Section b we consider some special cases, of the form *■ = v n 5 ,, , important 

lot the study of both "classical" and "exotic" hexagonal-cellular solutions (see 
Remark 6.1J, 


Since the vectors k in (2.11) are constrained by the requirement that 

f U * *S 

e have periods 2v/<t^ in x and 2ir/a 2 in y, it follows that n = I;.'., 

and m = k 9 /M 2 are integers and k must have the form k = (no, ,m - 9 ,0) . Thus, 
the only wave numbers, 0, corresponding to eigenfunctions having the required 
periods are those for which 


( 2 . 14 ) 


2 2 2 0 
a = n . rv“ + m utj 


9 2 2 9 9 

for some integers m, n, i.e., such that the ellipse a' = x ^ + y <:£ passes 
through at least one lattice point (n,m) 4 (0,0). (Note that there is no non- 
trivial solution of (2.32) if o = 0.) There are countably many such wave numbers 

U.' 'j >•) '• each of which corresponds to a finite, even number of lattice 

points ( ( n , tm) . 





1 


W't-j'ti » ¥ u 


x. *' »y 


H V 


l.< |( » on i«ti|Hm>ln to t lu* 2h ^ lot t ! re point;; 


(n . ,m ,) 
P] p! 


v - 1 " v l* "• 


v.*li‘* rt* tin* Indices an* chosen so that n , « -n ,, m , . 

p t- j ) Pi p(-i) 


- c Vi' llu n '■•'■ 


(.Mi) k 

an I observe that 


p| * Cn pi , 'T ,n p.1 i 2» 0) ' j * tl ‘ + 2 ”-» 4 V p * 1 ’ 2 ’***’ 


( l<>) k , - -k . . 

■i>M) p.J 

,1 ’ 1 each - - «» the reduced eigenvalue problem (2.12) has an infinite sequence 
oi real, nontrivial solutions 

! ‘I 3 f l ■ 1 1 » ' 2 , , . . . 

Hineo (- 1 tt‘y-4^) satisfies (2.12) whenever satisfies (2.12), 

wav order the indices so that 


C’. Id) 
The m 


p(-q> 


.M _ ,.pq .,,p(-q) e „,pq , , . 

pc| ’ "3 '3 ’ t\ h ’ 0 u P i p2 ' * 


pi 


are simple ei ge.nva 1 ties and the corresponding >,' p ^ ' may be cal en t* 1 


be positive on 


(-3/2, 1/2), Moreover, since, rn in (2.13) 


(s e(|iia 1 to in (2. 14) for sonic unique, positive 1 nt’oger p^, 


(2 , hO 


,: 1 “ V 


nun ]i 


Pi 


0 p=l , 2 , . . . 

Is also a simple eigenvalue of (2.12) and, for q )* 1, p > n^ if p f p Q , 
One now sees from (2.1.1) that the full eigenvalue problem / 2.10) has the 
so I at ions 


tor p <= 1,2,3,..*, q = ±1,±2 ,±3, . . . , where 


(2.20) \ * ]lp C] , v - i|- pql (x) = *P«( 8 ), ;j - +l,d:2, 


,f. pq - j (,) - 


1 <1 , , lit- , 

n ~ *pq -J. m A. A pq A pq , fi pq 

n 2 p;i da U ’ a 2 p.j da S ' ! 3 , ! / ( 

I 1 ‘ P 


(2.21) 
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to, 


\ou- that the depend on ,j only in the first erriponents. A* ■ - > r> i i : 

lu f ‘.in), (2.20), (2.21) and the fuel that v!^ ,f|'^ are real, w*. have 

i .•..*•) v |K|(H) - V" l] . 

it 1 m shown in the Appendix that the eigenfunctions ! ' in (2.20) • 

he assumed orthonormal in H , after scaling with constants depend inn or. p »nd 
q hut not on 1 . Thus we suppose that 


M) 


..pql rsu 

9 ’*' * 


A i* . 

pr qs Jt 


i.ht re is the usual Kronecker delta symbol. 

It 

The next lemma summarizes some of the properties lust discussed. The 
cnmi>, iciness properties are essentially well known (e.g., see [11]), while 
(2,,’T) is easily derived from (2,7). 


Lemma, 2.1. (i) The operator L :// ~ r H is hounded, linear, sel fndjoint 

.uni compact . Its characteristic values and eigenfunctions are given bv (2.20) 
and satisfy (2.22) and (2.23). The eigenfunctions are complete in H. 

(ft) The operator M: // ->- H is bounded, linear and compact. Its ad-joint, 
M" , is characterized l)V 

(M V 'v,w) » 'll Z'ew, , 


C 2 . .‘a) 


V,W € H, 
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'* * * '* 

5 * Tin* selection equal ions , We show next th.it tin* i’<*nt*r.t !!?(-*} i'en.rd 
iit.'Mi'i" i*. m ht reduced to a •' ini it— (1 ineiiHioual one. ibis r» duet i>*a in srci* 
it! !>v meins til spill lint* methods u*i I n?» t ho "si met ;n e" p.m limpet* . It o 1 . * > 
■ • tit ".it pJ it mli*" parameter. 

■ k '< ituv p. p, I veil bv (2,19) i s a simple eigenvalue ol (.*,12) .aid , / 

i ip 

• >r (p,q) 4 (p (i ,l), it la also a characteristic value ol I. of i-tult ii>ll» it v 

2N .’a . Tin* associated null space /( of 1 ■ t».I, is spanned t*v 

l ‘f) 1 

M.n J v p 0 lj , j - U,12,..., + N. 


< Ail. u deal {tip with quantities on il it will often be convenient to suppress 

Pn 1 5 i 

t If indices p “ p , q 1. Tims we write k . * k., i ~ ‘ , at* ,) 

J Pq 1 ' 

'IV » rt lui: «>nn 1 eomplomi at , J! 1 , ol M in H in spanned by ' : r P , q) , 

v/i- ilia 1 1 look for solutions of (*) having the lora v •* . + • with In 
M .aul in £1 , Tn order to studv the wav L. and F act on v it will 

I ! 

!m liieiul to inlrocluce some related operators. Let !':|f M 1 denote the or- 
tic .mal projection ol H onto M and let KiM M‘ denote the inverse oi 
the restriction of I - p^L to M 1 . In addition, we define bilinear operators 
:li II * II and J‘:f/ * II -> 11 by 

u,v,w It 
u,v,w H, 


l 1 . * 


c 5. 1) 


1 ; (tl ,v) ,w) = - (u*7v) • w, 


wv 


{i(ll,v),w) “ 

Om* sees easily from (2.8) and (2.9) that 


(M) F(v) = T(v,v) and 0(v) = f(v,v), v . H. 


It will frequently be convenient to represent v t li by its Fourier series 


( i. r i) 


v 



pq.i 


j 


where the sum is extended over the set of integer triples (p,q,j) with 
! ' jj • N, 1 1 p ”• . I '"|q| * . Wlien v '■ ; . M, ( J . r 3 ) becomes 
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( i.Hjl 


I 

h’-i 


-mi 5 _ i»m (« I ) 


It follows that v In il.'i) (.ir In < I , n 0 • 1 , < 


>ll * W-1> ( " r ‘l * -.I 1 ' 


■ » f l! t.' 

It mu nil 1 v ii 

Tin* J ol lowing lemma, proved In the Amend lx, tumbles us in Horn* .situations 
ii> . ilmlate with the operators introduced above, Here and throughout the pa. -a 
( 5 is zero whenever the (scalar or vector) parameter f is not zero, ami to) 


I ci’inui. 1. 1 . (1) li H, then Lv, Mv and K Pv can be obtained !Mt 

l>v formal calculation, E.g., 


(5, ) 


1 v » , KPv * f .-m- • nni 

IWIVIW ’ ' ,r pq.1 l> pq -»y 


e liei * • () denotes suinmat Ion over me same set of integer triples as ! , e:-.ce, t 

that (p,‘i,j) ¥ (p () ,l,.j). In particular, K is bounded, positive and ■>. 1 fau'-nint 

mi 

H !' i'JSM * '1 1 , l.e., for < M, * 0, jij ~ 1,2,,.,, N. 

UMJ For u,v,w ■ II, (i(u,v) ,w) « ~('!(u,w) ,v) and I’(u,v) * ”(v,u), 

( iv) 1:11 ' M > .M x . Tn particular, F:M < M 1 and ( . (v ,lr!1 , . ' ) ," n ) a 0 

fot ill p,q,r, l,n with 1 ■' p 1 * j q J : '»* and | r| , t i ‘ * !n* * 1,2 N . 

tv) If has the form (3.6), then there are real constants b fj ,h,,b,,bj 
depending on but not on n, such that 0, 


( }.«) 

( i . y ) 
0 . 10 ) 
and 
(3 , II) 


(MKMz,T n ) - b (J B_ n , 


(MKK(ilO ,7”) » b 


=1 ;VV 5 % + k r + k n> * 


(*I'(i!',KM ! ) , 7 ") - b 


(oOlO ,7 n ) » b. 


I * r 


N 

?■ 6 , ?> ‘ (k. + k + k ) 

,|r|=l J r -- 1 r " 


,W% + k r + V- 



I 5, 


(vi) There are nonnegative constants a . such that 

P/j 1 11 


( 5 . ;.») 


(;•( ,Kr(, )),~ n ) 


j: a , (2 - \ ).-,;• .e 
I l-i P 0 i° i» 1 “ n 


t hs> const nit 0 ^ j () depends only on Pq and * 

,i . " a , * a. k , w The exceptional eases in which the constant 

P {) (**nK~j) 

i i a ’em are described In (A. 24) and (A. 17) of the Appendix, (Met- ,il .o 

ii ii 11 

’ o * 

?:• ■ ; ft ■ . I 1 . 


Ikj + k n i so that 


Uv .hall need to relate the spectral analysis of the linear operator 1, 
to that of the linear operator L " I. - >M. For small values of . it is well 
In own (e.; r ,,, see f8]) that the characteristic values of 1^, are perturbations 
ni i hose o! (, , In fact , the characteristic values el L,. are determined hv 

i 

tl’.e prolilem obtained from (2,12) upon replacing (2.12a) with 

(1.11) 0 - 1)S 3 - .Vi 2 (1 - 2 -/k)=J» 4 . 

Dm* ‘ (mis, in particular, that the critical characteristic value, ^ * '^(y), 

i.e., the characteristic value of L, of least magnitude, is real and simple 

V 

as .a eigenvalue of tne problem (2.12) with (2.12n) replaced hv (3.11) and 
.cl c.-jtia1 to s () . (The relationship between and the critical Rnyleldi number 

is tin* usual one described in ri;4;7].) The next lemma specifies the expansion 
in . of \ (i and niav be proved along the lines of the development far the non- 
I fin 1 . n problem leading to equations (3.17) and (3.18). 


l.enmia 3.2. The critical characteristic value, > , of L., has the expansion 

* I 

(3. i'i) > c « Uj - r' 2 lijb 0 + A c (y) 

where p. is given by (2.19), is as in (2.8) and A^(y) is real and 

| *‘ 0 (Y) | » 0(y 3 ) as Y •> 0. 


su t i sl ies 
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I'ur small i we seek u solution of equation (*) in the form 

N, !>) v * ,( + Vi). >■ * ** - ). 

Sit ri- v .M, T * .‘f 1 and t < R^ art to be determined, T f « () + Of.), 

i lien x - ** * + ) and a .solution of tlu* form (’3.15), for small 

ia ailierll leal If 0 or superrrit ienl if *- <>, 

We aiib.sLlt tile (3,1)) in equation (*) , use P and 8 « J « P to pro jeel 
.into M 1 and and use (it) and (iv) of Lemma ? to obtain the following 

equal i euti on M* and M: 

(*."»> (.a) 1) - (i - UjL) I + U r M i- PC) + i P [ , : } M * - C.VO - ’C, ■) - G(.j] 

+ T ? F , f-n ] (uJb 0 - r)M‘," + H^l^bo - UL ! - V( ) - 2’( > )] 

+ , 1 P[-n 1 (i'*jb 0 - ' )M‘:' - «('?)), 

(b) n « C' 2 i> () - <),- + sliijM'1 - •!■(', )-{;(;]+ YSfp(0 - ;i (.,)! 

+ i' 2 S f-)J , (jj jb Q - i )M ,; - R( •)]. 

Since K -■ (I ~ u^l,) * is bounded on M 1 , given 1:^ *• 0 tliere is a 
r 0 such that if (,,i) * M * R 1 with Jrj + j| ■:•}! »■' t () then one ran selvt 

( 5 , Ilia), by successive approximations, for ’I' ~ whenever **• j *' < „ . 

in ! art satisfies 

(3. l/> Y * -p j KM 1 }' + KI’C'O + Y'f-, , 

where Tj « i^C,, <',i) • M 1 is bounded depending only on t^. We next use (3.17) 
to eliminate V from (3.16b), taking (3.8) into account Lo get 

O. 18) 0 = -T. + S[p 1 MKF(t) + U^'O'sKM*) - (<!•«)) - G(v)J + RC , , ,') • 

Here, for jy| < Yq and |t| + |j^|| < tg, the remainder term 

ROfvivy) = ys (i* 1 M'l'j - ‘K'Mj) + V(V - 2P + vHi/nfbo - 


€3. 1*1) 


) M ,; - (’<( )) 
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satisfies, fur Home ■ 0 depending only on , 

( 1 • 2U) 1 1 R (>' » l * k ) ij ^ V £q • 

We take the inner product of ty 11 witli equation (3.18) » making use of the 
expansion (3.6) and various formulas in Lemma 3.1 to obtain 


( j.2i j 


0 = h ( 3 * i » r) ~ -10 , + b 


-n 


! 1 i . : J ! “i 


k , + k . + k ) 
1 j i j n 


+ I j - i 1 1 > 0* n ^ 2 ~ \jn )k V''-.i 1 '- n + 'n 1 ’ r) ’ l’*' " ‘ 


Hete b = li j(bj + b 0 ) - and, according to (3.20), 
(3.22 r n ( n,i,i) = (R(W, t ,'i) ,f ") , 


sa L is I Les j r (!•» > ,y) ! YFq - 

For the reasons discussed in the introductxon (see also the discussion 
in [ P j ) we must augment the system (3.21) by an equation, ('(r, O ~ 
Involving the so-called generalised dissipation V , where r is a 


real parameter and 
(3.23) (/(.;» 0 


1 N 1 
4 I • lijrl , + lb 


1.11=1 J ' J 
N 


3 • . 


i i . 1 i , ,TJ m i 
! » I j ! » ! m i ==1 


v- . ;* . .* ’ ( k . + k . + k ) 
’ ■ j m 


X —» V a (* 9 — \ 2 (»* O J 

4 ji| , | j |=i Po iJ iJ ‘“J-i 

rims, we consider the system of selection equations 
(3.2',) (a) 0 = F(,’>, i ,y) , 

2N 1 

(1>) - I'd:, i), r- • t , (r.-y,,) c IR , 


'** n f " (f n-’W-l N and » * < S -H B -Pl V- 
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Tile functional \J la essentially the functional K in [3, p.631] 
with t - y. In fact, setting e ■ y in the analysis in [1], one obtains 
formally a number of expansions, equations, etc., that are closely related 
to various quantities used in the analysis here. 

Perhaps one would hope to solve (3.24) by solving the equations, e.g, , when 
(Y,< ) - (0,0) and then using the implicit function theorem to extend such 
a solution to a small. (Y,t') neighborhood of (0,0). One anticipates, however, 
difficulty here in implementing the Implicit-function theorem argument 
(e.g.,, see [17]) because the equations are invariant under translations 
of the (x,v)~plnne. Consequently, the solutions will not be Isolated 
•iml the relevant Jaeobians will be zero. Thus, it is natural to seek 
solutions in a subspace of H, where one may hope that solutions will 
be isolated. This is conveniently done in the next section in terms of 
group representations as in [17]. 
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4. I'Ju* reduced sel ucLlnn equations. XI k* bn sit* subspnoo , $ * used 

throughout the ruma Inder of Liu.* paper is introduced in this section to- 
gether will) some technical lemmas regarding real solutions of equation (*) , 
hut v be the 2*2 matrix of a plane rotation or reflection and 1 « t 
a « be a translation vector. For k « 3,4 let denote the 

k * k matrix ob.Lai.ned from the identity by inserting r in place of the 
2 2 Identity matrix in the upper left-hand corner, Set * (aj,a,,Uj 

and let o « fr^.a^} represent an arbitrary plane rigid motion of x - (’•;,y l z) 
space that keeps x fixed: -x * r~x + a,,. then a representation, • - 1, 

oi this group, (i, of rigid motions is defined by 

14.1; (T,v)(x) 3 r^v(o *x) 

lor smooth tour-dimensional vector fields v defined fur x . Bl^. 

When y « 0 it is well known (e.g., see [11, 17]) that the Boussint-sq 
equations in (2, 1) are invariant under T t for o ■ G. The next lemma 
shows that a corresponding invariance property holds for equation ( * 3 when ■ = 
and that Lhe invariance also extends to the case ) 4- 0. Such an in- 
variance statement makes sense, of course, only for eg v for which 
both v and T v J ie in H. 

' 5 

l.e nnna 4.1 . but o G and suppose that u, v, T.u, T,v all lie in 
II. Then each of the operators L, M, i, I is invariant in the sense that 
l,(T..v; = T } (hv), 1(T u,T,v) = T i’(u,v), etc. Consequently, 

(4.2) L (T v) = T (L v), F (Tv) = T,(F. (v)) 

\ / y a o y y o o y 

so that equation (*) is invariant under T . 

Proof . Kach of the operators L, M, I' is defined ( (2. 6) , (2 . 7) , ( 3 . 2) , 

/ 

1 1 . 


(3.3)) by an integral of the form 


Aw, where A is a linear, A(v), or 
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l [linear, A(u,v), term in the Bnussinesq equations. Il A is invariant uml> r 

T . , i In n 1 t Is easy to see that the corresponding operator is invariant under 

T . li . r . , lor as defined in (3.3) note that A(u,v) (OjOjU/V, ,0) , 

a a 

Since (r^).,j * ‘j.j, A(u,v) is invariant under T, because 

( i ,A(u, v ) ) f x) “ (0,0,u^(o ^xJv^O - . ^x) ,0) * (0 ,0 , (i ,u) (x) (1 v),fx),0). 

The invariance of Au « (0 *0 ,2xu^ ,0) , corresponding, to the operator V, ■'el1ov> : 
in a similar way; the invariance of the A's corresponding, to 1, and l« 
proved in f c ) j . 

Remark _4.1 . Because of Lemma 4.1, we may study problem (*) on any of the 
rinsed linear subspaces, S - \v r l(:T,v - v), without the use of prelections, 
bv Merely restricting t’’o operators in (*) to S,. Under such a restriction 
the equation (*) is denoted (*) and retains its form; similarly the new sel.-e- 
t i lit equal ions , <’3.24)„, are obtained from (3.24) merely by restrict inc the uo- 
effteienis r in a well-defined manner determined by T . By the restriction 
to $ we shall avoid the problem of zero Jacobians mentioned above. Of course 
n solution of (*/ (T in is also a solution of (*) in H. On tin- other hand. 

;i Kiiibilitv proof in S n , although encouraging, is a weaker statement tlian one 
on II but inst. iluljty in 5, does imply Instability in II. 

Throughout the remainder of the paper we shall largely restrict our atten- 
t ion to and its subspaces, where S n • when 0 denotes rotation by 

n radians about the x-axis. Thus, (x,v) -> (-x,-v) under tt and 

(4,3) (T. [t v) (x,y,z) = (- v 1 ,-v 2 , v 3 .v^) (-x,-y ,z) . 

Tt follows fron (2 . 20)- (2. 22) and (4.3) that 

T (i p<iJ = ,jji ,c l <" i ) = . 


(4.4) 
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I v , v « '%(’ in (3. 6) 1 1 or In ,S' !t it and only if e.n h m the 

i’ih* ill eien t s sal Is! ins i ‘ Jn fallows direct !v fre* 

(4.4) that T sat Is) ies 

(4, >) (T, t u,v) * (u,T 7J v), u,v < H, 

Will U> II is a complex Hilbert space, we are of course interested onlv In 

real .solutions of (*). Since the basis elements satisfy (2.22), the om*f t l*-i, rut s 

In i he expansion ( 3 . r j ) satisfy F ^ = i* , whenever v is real. In i ..ri loi- 

l.i! , !->t real v In -S equation (4. r )) Implies further that the . . ar* 

P'l ) 

ii'i 1 , Moreover, the following, lemma shows that the operators !.• equation (*i n<* 
real operators. 

la ■iniiui 4.2. (i) The operators L, M, 1, F are real in the sense tlr t 

l.v l.v, ;*(u,v) = ’Ku,v), etc. 

(ti) If ! > < M is real and y,'t • K* satisfy |vj ■' Tq, : t | + ’’ t () , 

with i ( j and t () such that (3. 1 7)- (') . 20) hold, then f - :(,;,.) in (3.17) 
mil It It in (3.19) are real. 

Proof. Part (I) follows easilv from the def In i t ions , since the ret respond I n 

d i i i iM'ential operators have real coefficients, e.g., (F(u,v),w) - ( r (u,v.) ,w) ~ 

* 

j u v /( w^ = (F(v:,v),w). For part (ii) , note that if y,T and if are real then 

upon taking the complex conjugate of (3.16a) and using part (i) we see that ; 
is a solution of (3.16a) whenever ¥ is a solution. But the successive-approxi- 
mations solution of (3.16a) is unique in a small neighborhood of + KT( ) , 

which is real. Hence ¥ = ¥ is real and by (3.17) ¥^ is real. From (1) 

and (3.19), R(¥,I,y) is real. 

Since, according to Lemma 4.2, ¥(i|t,T,y) is real whenever y,T and 

>!’ ■ M are real, the problem of finding real solutions of (*) is reduced to that 

2 

of finding, for sufficiently small (Y,£) c R , solutions (S,x) of the 



selection equations (3.2A) with 

<A.<>) ! < Hi 1 , - P.,, 
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N i 

t and ■ J! p iH 

h!-i J 

.1 ** 3 ,2,. . . ,N. 


real, I.e., 


wi th 


20 , 


In the remainder of this section we consider problem (*)_ obtained by 
restricting (*) to S, n . The nullspace of I - v^L restricted to S r is 
‘I,. -M n S. { . From (A. *3) we see that if t{» <. then ^ « p . , j * 1,.,.,N and 

N , 1 N 1 

(A. 7) T * /. ti . * K P . (p + >l> • ) 

! ,i |*i J j-i ^ 

Thun, M. is N-dimensional and we shall henceforth take the liberty of sup- 
pressing »’•_ k_nj in the notation, i.e. we write • = instead 

oi - .P- .... ,P M > and we regard f and V as functions of 

n j i in n 

N 2 

( , ,[) tn (E Moreover, in the context of S we have the following 

lemma (see also the related results jv) [17]). 


1 

Lejimiu _A_.J3_. If ’/< »: M and *y,.T •' E are suf ficiently small then 
[ ( ,i,r) F , r ,y) , n = 1,2,. ..,N. If, in addition, I is real then 
r n r n (i,!,y) in (3.21) is real, n = J ,2, , . , ,N, 


Prop / . Since 


N 

= B. whenever ifj ~ £ B.,^ belongs tn )(„, , and 


since a . - a , w .. in (3.21), to show that F - F it suffices to 

P () nj P 0 (-n)(-;j) n -n 

show that r = r in (3.22). Using the fact that i\.\‘ = for ‘ . M„ , 

one sees from the invariance of (3.16a) under T and the uniqueness of "t that 

T - T also holds. It follows that T H'« = H*, and T R = H, where *. and 

t ti 1 1 T 1 

1! are given by (3.37) and (3.19). Thus, one sees from (A. A) and (A. 5) that 
r - (R,!' n ) = (R.T^tj; n ) = (R,i|j n ) = r^. If, in addition, ^ is real then (2,22) 

imply = (R,ip” n ) = (R,tp n ) = r so that r n = 


and Lemma A, 2 


r . 
n 
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Because of Lemma 4.3, the selection equations 0.24) in the sett in/ of 
M mav be replaced by an equivalent system of N + 1 equations In the N 
(possibly complex) variables b * (bp...,b^), the real variable i , and the 
real parameters f and / : 


W) o - f ovi.y) -in + b z A t1n iy , 

1»J“1 


+ }. A, i '•> + r (b,t,f), n * 1.2,-, ..,N 
. 1 in i n u ’ ’ 
i*l 


00 


N 

t'OvO •; -i h: + 15 


N 


A, 


1 


N 


o > 


. ■ i ■- ... , [‘i.l . r< *1* ?, A , , > . ; . , 

) 1 . . , t in i m 2 . . 1 i ] i 


wlieri* (b, 1 *'?'»*■-) * , 


(/,.<)) ft - *(k 1 + kj + k.) + ^K+kj-k.) + Mkj-kj + k.) + *<k, - k, - * 


0.10) 


A. , = a , .(2 - 6..) + 2a A / A \ 
ij P Q ij P 0 i(“j) • 


Moreover, since Lemma 4.3 shows also that F = (Fp...,F^) mav be regarded 

N 2 n 

as a mapping of a neighborhood of (0,0,0) in E v E into E , it is 

natural, to seek solutions of the selection equations in (4.8) of the form 

N+1 

(:*(.,0, r* (y,t) ) > IR by use of the implicit function theorem near 
Y - 0. Tf (/*,i*) < 11 N+ ^ is such a solution of (4.8) near f = • = 0, 

then (b* . . . ,b*,b*, . . . ,I3J,T*) is a solution of (3.24) satisfying (4.6) with 
; . - [’ , , i.e., a solution of (3.24)^ satisfying (4.6). Thus, the above 

construction leads to real t|j in M and, hence, real solutions of (*)„. in 
,S’ . To actually carry out the above construction, we seek first the real 

‘H 

solutions of the reduced selection equations obtained by setting Y - »*■' = 0 
in (4.8): 

(4.11) (a) 0 « F n ((i,i ,0), n = 1,2 N 

(0,T) c E N+1 . 


(b) 0=l/(fi,T), 
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Remar k 4 ♦ 2 , It is easy to check that ( c' , r ,0) * - ’iiifj'il t „ « \ , 2 t . 

so that (4,11a) Is a gradient system. Since (4,11a) Is not the reduced bi fin 

cation system associated with (*) , this gradient structure is not identical 

to that used extensively in (1, 10, 173, although it is closely related, We 

note that the reduced system obtained from (3.24) by setting r * , * 0 has 

a similar structure, with F (p.,',0) * 1/(15,*')} the factor 0 appears in 

n 

Llie 5, case because of the identification of and . 


In developing a selection principle for stable suborltical hexagonal 
cells one needs to consider only the. reduced selection equations in (4.11). 
Other choices of the reduced selection equations are also appropriate in 
convection problems, e.g., in the study of supercritical solutions and the 
exchange of stability between rolls and hexagonal cells, and will be con- 
Hidered in a subsequent paper. 


r* / 


ORKHX/H, t ..r ■ . 4 
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5. KxlBti‘mn‘, and Htability of rtMlsolut ions in .S . Tn thin sect inn 
we solve the selection equations In a general setting bv moans of variational 
methods. 

The t allowing preliminary result yields real solutions of (*) in 

Theorem r \. 1 . bet (i *,''*) • 11 , with * y p, satisfy the reduced 

select ion equations (4. 11) and suppose the Jacobian det jeH 
(: , ,1,0 « (!'*, i * ,0,0) . Then there Is a o 0 such that for (,, ) R‘ 


is not acre at 

J 


with 


<S the selection equations (4.8) have a solution 


N+X 


(»■ (-,0,i (,, )) IR satisfying 


0 . 1 ) 


lim (t’(Y,0 , f <'y **■)) * (••*,!*) • 
(h e ) on,o) 


i‘u rt he rmore problem (*) has a real solution of the form (3.15) wi th 


N 




" (yv ) (T 1 + <{> ] ) » > ** r (y *f ) » 


>1 


and '1' obtained from »j»,T by means of (3.J7). 

The result follows from the implicit-function theorem applied to F,(/ 
near (;’ , i ,Y,t ) * ((•’*, r*,0 ,0) , provided that dot is not zero when 

evaluated at (ft*, i *,0,0) . But “ 2 F is zero at this point and 

tip 1 i 

Thus , 


& i ' 


i- 0 , 


dot (P*,T*,0,0) - - 1 8* j 1 det 


4J 7 ] (i’*,t*,0,0) i 0 


!)({■!, t) 

and the rest of the theorem follows easily. 

To utilize Theorem 5.1 we seek solutions of the reduced selection equations 
with f‘>* i- 0. We next show how this may be accomplished by exploiting the varia- 
tional. structure of the reduced problem (4.11). 


Note that 
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0.3) 
where 
0.4) (a) 


Jl'GvO - - i r|i j 2 + q(3) + f(fc), 

*» A 


qGO 


(b) c(i») 


3 i.lVm-l n ' 

, N 2 2 

4 V A p. (* 


In order Lo determine suberltlcal solutions of (*) we shall impose the fnllnwi!« * 
hypotheses on q and c : 


on q(}‘) / 0 on 311^ 

(IM c(, ) :■ 0 for all I) 4 0 in K N , 


RomnxliJ_.,l ; , Hypothesis (H^) fails, in general, since -Kk^ + k^ + is 

zero unless the vectors k, , k . , k form an equilateral triangle: k. 4*k, +k * 0. 

"'i -.1 m -.I j m 

This latter condition is possible for hexagonal lattices, oj. * /3\, > n * /, 
when i> satisfies (2.13) for integers n^, m^ of the same parity. In such 


cases (ll ) is sat is fled if b 4 0. Concerning (H e ) , the condition j “• 0 

follows from (4.10) and the nonnegativity of the a^ ^ in (vi) of Lemma 3.1. Ho 

hypothesis (II ) is satisfied, e.g., if a . >0 for^each 1 * 1,2,,..,N. The 
(. \y () xx 

latter condition is fulfilled if at least one term in the sum defining a 

IV i 

is different from zero, 

In the following discussion of the finite-dimensic-nal problem (4.11), a 


prime denotes the gradient with respect to ft. Thus 


c 


! 


(:•) 



N 

. 1*1 


f"(B) 



,J-1 


etc. 


In view of Remark 4.2, the system (4.11) becomes 

(5.3) (a) 0 = -X?, + q’(6) + c’(3) , 

(b) 0 = - £|B| 2 + q(B) + c(B) . 
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Wo deli no selection, t uncjJjnyU s f and g by 


( f M>) (a) 


(aOa + rCO _ lf ,, 4 o 
j ft 

! o , if t - o 


* 


00 K(* ) «| 2 (| )/(4e(l )) . 

LiMiuta^ 'i , l » l.et tin* functionals q and <■ be given by (570 and supper-.- 
i' sat Ini i i*B (H t ). Let G‘,t) * 1 with i t Q and net f d ,7 j ! . Than 

tlu' lol loving are equivalent. 

(i) (3,0 in a solution of (5,5), 

•f 

(li) B is a critical point of f with critical value f(h) ■ g , 

(ill) B is a critical point of g on jdj * 1 with critical value 
}*.(■') ~ ** , and tlu* magnitude of !’• satisfies 


(5.7) 


hc 1 00/2c(li)l. 


Proof. The critical points of f(E) are determined by 


(5.8) 


0 * f'G‘) * + q'CiO + C* (B) } , 


where 

(5.9) -t * 2 [q (I ) + c ((•;)] | H |”“ 2 * 2f ((•;) . 


Since B i- 0, equations (5.8), (5.9) are just (5.5). Thus (i) and (li) are 

A I I 

equivalent. The condition that B be a critical point of g(6) on jiM * 1 

T 

with critical value - 2 is 

(5.10) -IB » g'(B) - q(3)[2c(B)]" 2 [2c(3)q'(B) - q(B)C(6)] . 

If we use the homogeneity of q, q', c, c’, g' and the Euler identities 
B»q ’ (ii) « 3q([i) , (>c , (B) 83 4c(B)» B*g'(3) = 2g(B) , then from (5.10) we get 



(5.11) 

( 5 . 12 ) 
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-i * i-.-g'O ) * 2 ^ (: ) , 

-if' * jvJg'CiO * ft* Ci > - q(l:‘) ) ]” 2 I2c(i0q' (»0 ** q(;)i‘ ! (. )). 


If (F,i) satisfies (5.5) with i i 0, the* upon multiplying (5.5a) by t, 
using t lu* Hu let* Identities find subtracting twice (5.5b) we obtain 


(5.11) 0 - q( ) + 2c(iO, 


which iiupUaH (5.7). From v 4 0, (5.13) and (H^) we have q(t ) « -2c ( ) 0. 

Fut . itch equations (5.12) and (5 . On) art* the same. Similarly, (5.5b) and 
(5,11) imply 

2 


‘ i* |" - [q (! ) + *'({'<) ] 


mo iliac (5,5b) and (5.11) are the some. Thus (i) Implies (iii). Finally, let 
\i and | i,l satisfy the conditions in (iii). Since ,■ 4 0 bv assumption, 
(5./) is equivalent to (5.13) so that again (5.5) is the same as (5.11), (5.12), 
Thus, (Hi) implies (1). 

It is clear from Lemma 5.1 that solutions (b*,T*) of the reduced selec- 
tion equations with E* 4 0 are obtained from those, critical points /* of 

A 

g(t) on \ti\ » 1 for which g(B*) 4 0. Fuichermore, it follows from 
(H ), (5.6a), (5.13) and (ii) of Lemma 5.1, that for such critical points 


(5.14) r* - 2F(H*) * -2 - c -^- < 0. 

|H*| Z 

Thus, on the basis of (3,15ff.), a solution of (*)., generated from (f*,T*) 
will be sub critical* at least for small values of y and r . According to 
Theorem 5.1, to extend such a solution of (5.5) to a solution of (4.8) we must 
show that det ^ 0 at (B,T,y,c) = (B* ,T*,0*Q) , i.e. det h 4 0, where E 
is the symmetric matrix 



('>.13) 


* » r • 

IWH4 
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.r«v 


K - -l*| + q"(; *) + e"(H*) . 

Thun, dot 4jj! is zero if and only if E is singular. We have* established 
the following result. 


Jb 2 • Suppose q and c satisfy hypotheses (H ) and (H ) . 
Let (t*,t*), H* i 0, be a solution of the reduced selection equations (5.5) 
Kiteh that the matrix E in (5,15) is nonsingular. Then there exist ^ 0 

and i\ > 0 such that, when J Y J < and |f. j ^ t ^ , equation (*) 1T has a 
real, subcritical solution (v*(Y,f ) ,>*(y,t)) of the form (3.15) with 
i « i(y,0 < 0 and generalized dissipation (/“**. In fact, 


(5.16) (n) v*(Y,f ) = Y " 3*0p j + >f j ) + V(y,0, 

j»l *' 

(b) X*(y.>) * - Y 2 h 1 (y|b 0 - 't*> + A(Y ,0 * * c + + '( iY ) 

2 o 

where i* satisfies (5.14) and, as y -»• 0, V(Y,k) = 0(Y ),..(Y,e) = o(Y“) » 
/'(,-, ) - o(Y 2 ). 

According to Theorem 5.2 and (i) and (i i) of Lemma 5.1 we can generate a 

N 

solution of (*), t by finding a global minimum of f on H‘ . If (■’ = 
with ji’| * 1, note that 


(5.17) 


f(3) = c(3) 


P + 


q(g) 

2c(3) 


- g(3) 


We minimize f(P) on E by choosing p = -q(3) /[2c(P) ] and maximizing 
g(P) on 1 3 1 = 1. If q(3) £ 0 then we generate in this way at least one 
nontrivial solution of (5.5), say (g*,T*), with T* satisfying (5.14). 

If we differentiate (5.8) and make use of (5.8) and (5.9), then we find that 


(3.18) F"(3*) » |3*r 2 |-T*I + q"(R*) + c"(3*) ] * 1 3* | 2 I2 . 

^ince f has a minimum at 3*» we know that f '($*), hence E, is at least 


positive semi-definite, 



(5. 19) 
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'i'huH , If ]■) is nons ingular nr n minimum of f, it must he* positive* definite; 
wt* shall see tlmt the solution of (*) 7f generated (ns in Theorem 5.2) from 
(i-'rt.irt) is then stable in 5 , 

7T 

The relationship of the critical points of f to the generalized dissipation 
I* is given in the following remark. 

Re mark 5 .2. From (5.3), (5.5) and Lemma 5.1 one sees that if is a 

critical point of f(h) and r 0 - 2f (ii Q ) , then R () is also a critical point of 
t / (iQ,l>) and ^ ( * q , i'q) * 0. If Tq is also the absolute minimum of 2l'(r) 
then 1/ » 0 is the absolute minimum of U(Tq,6). 

We turn, then, to the question of stability of a solution (v* ,.'*) * 

(v*(f,0 , >A(V»^)) of (*) n having the form (5.16). For small y and ", the de- 

rived operator, V , of (*) at (v*,A*) is a linear Fredholm operator of index 
zero, the perturbation by a small bounded linear operator of the self-adjoint 
operator T - ]j^ L . As observed in [17], because of the invariance of the equations 
under translations of Lhe (x,y)-p.lane, the stability of solutions of (*) in 1/ is 
always indeterminate, In the case of S 77 , however, we have the following result. 
(The notion of stability here is "linearized stability" as in [ 16 ; 1 7 ] . ) 

T heorem 5.3 . For y,e sufficiently small, a solution v(y>€) of (*)„ obtaine 
from Theorem 5.2 is stable in at \ - X ( y ,c) if all eigenvalues of the matrix 

12 in (5.15) are positive, and unstable if some eigenvalue of E is negative. In 
particular, if v(y,F:) is generated from (3*,!*) corresponding to a minimum o„ f 
and such that E is nonsingular, then v(y,e) is stable in 5^ at A = A(y,c). 

To prove Theorem 5.3, one proceeds as in [16] to determine a subspace of 

5 invariant under V and corresponding to the N critical eigenvalues of 
V for sufficiently small Y and £. This subspace has a basis of the form 
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+ yz 


h: 


1,2 N 


sat Is Tying 

1 ^2 1 

( 3 . 21 ) Vz l » /' Y b..z , 1 - 1,2»...,N. 

. 1-1 •' 

To establish the existence of the basis (z , ...,z N } in (5.20) one needs to 
show, in particular, that T. Z* * Z* so that ?} belongs to 5. . The proof 
that T Z* * '/} makes use of the fact that (ty* + ;~ i ) belongs to ,M . and 
ful lows along the lines of the derivation of (‘3.17) and the proof of I.emmn 1.3, 
hinee we assume in Theorem 3.2 that K is nonsingular, (5.19) implies that all 

eigenvalues of E are positive if E* minimizes f. Tn this case v(j, ) 
is stable in S,, at > * >. ("w) for (y,r) sufficiently small. 
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6 . Subcritlca l hexagona l cellular solutions . We now restrict the 
problem to the hexagonal lattice and prove a general result about stable 
‘iuberllic.il solutions which yields a selection principle for hexagonal 
cellular solutions. To fix the ideas we treat also the special case 
of dim M = 12 in Remarks 6.1, 6.2 and 6.4; this case is the 
selling in which "exotic" solutions of the Bennrd problem were originally 
studied in [ 10 ]. 

We begin by showing, for an unbounded sequence of Integers N, 

that one can determine N sextuples of critical wave vectors corresponding 

to the critical wave number, Uq . These 6 N vectors generate a nullspace, 

M, with dim M = 6 N. Take = /3a, ~ l and choose Ct so that 

N 

(2.13) has exactly N distinct solutions (n,,,m n ) = (n.,m.). where 

0 ’ 0 j ] ]=1 

n. and m. are nonnegative integers of like parity for which the critical 

wave vector k. “ «x(/3nj ,m^ , 0 ) makes an angle h 0 < '•^<””/3, with 

( 1 , 0 , 0 ). I.e. , take a = Oq/v^Fq where the integer Mq is chosen so that 
2 2 

the equation 3n + m = has exactly N solutions satisfying the 

above conditions. (Tt is well-known that such pairs (N,M^) exist for 

an unbounded sequence of integers N (e.g., see [20, p . 345 ,ex. 5 ] ) .) 

We suppose the N vectors, k^ , are ordered so that 0 £ 0^ < 62 < ... ^ 0^ 

Define the N triples, T_. = (k^ ,k^ + ^ where, for j - 1,2,...,N, 

k. +N (resp., ^ 4 . 2 ^^ obtained by rotating kj counterclockwise 

through rr / 3 (resp., 2 tt/ 3) radians. Note that the 3N vectors, k^ , 

have lengths and direction angles Ch satisfying 0 < 0^ < (.^ < * • • < r ' 

Bach of the N triples, T. , can now be extended to a sextuple, (T, , - T. ) 

J J >1 

if we define k_^ = -k^ (j = 1,...,3N) in accordance with (2,16). In the 
above context there are . '^finitely many possible period rectangles 


corresponding to values of a = however, the critical wave number, 


31 . 

Op, and the "size" of the basic hexagonal cell remain fixed throughout 
the following discussion. 

R emark 6.1 . If In the above (N,Mq)*= (1,6) then n Q « 1 and 
nip * 1 in (2.13) and dim M * 6. In this case we have one triple 
(k.| ,k,, ,k.j) and one sextuple Ckj »~^2 * w * ier ‘ a ^ = <>(t r 3,],0), 

k 2 * <(0,2,0), and k 3 - <t(-/3,l,0) f If (N,M Q )= (2,28) then 
lip = 3 and mp = 1 in (2.13) and dim M « 12. In this case we have 
two triples (k^.kg.k^) and (k^.k^.k^), where 

k- « c*(3/3,l,0) k 0 - a(2/3,4,0) 

(6.1) k 3 « a(/3,5,0) k 4 = a(-/3,5,0) 

k r = a(-2/3,4,0) k, = aC-3^,1,0) . 

J ~ o 

The first of these special cases, dim M = 6 , was studied in [ 2 ; 5; 9; 17] 

In the context of classical hexagonal solutions. The second case, 

dim M = 12, was studied in [10] in the context of "exotic" solutions. 

6N 

We now define a basis {( b.}, n for M in accordance with (3.1) 

j j = i 

and proceed as in Sections 3 through 5. To make use of Theorem 5.2 in 

3N 

the present setting, one needs to minimize f on B , where f is 
defined as in (5.6). Thus, we require, in particular, the coefficients 
in the functionals q and c defined by (5.4) with N replaced by 3N. 

The coefficients of q are given by 

1, if (i,j,m) is a permutation of (n,n+N, n+2N) for 
some n e (1,2,... ,N} 

0, otherwise. 

Thus, setting B = 2b with b defined as in (3.21), we find 


( 6 . 2 ) 


A. . = 

13 m 



r» r 

'U . 
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N 


(6.3) 


t](p,) * B ,^ ^j+N^i+atr 
3 ~ J 


Note that if b ^ 0 then q(3) £ 0 since, e.g., - ^j+N + ^ j +2N * ^ 

lor j » 1 ,2 , . . . ,N. 

Wo next discuss the coefficients a^ and required to determine 

o(i) (see (4. 10), (5. 4) and (A. 37)). For our purposes it suffices to 

evaluate a., and a., when k. and k. lie in the same triple T , 
ij i(-j) "i *»j 1 n’ 

n * 1,2,...,N. Recall that depends only on jk^ + k . j , i.e., 

onlv on the angle between k, and k. (see (A. 37) & ff. ). Wien k. ,k, 

1 •’0 -.1 

lie in the same triple and k. 4 k. this angle is either tt/ 3 or 2~/3. 

»v X *v j 

We denote the corresponding values of a^. by a("r/3) and a(2~/3), 

respectively. It is now easily seen that if = a (tt / 3) then a^_.p = a (2 "/ 3) 

and if i.. = a(27r/3) then a., ,v = a(ir/3). Since A.. = 2(a.. + a , , .-.) when 
ij i(-j) ij ij i(-j) 

i 4 ij , and since a , ■> = a . ✓ . \ when a = a.., it follows that 

p(-q) i(-j ) pq ij’ 

the A„ have a common value, A = 2(a(ir/3) + a (2^/3)), when i 4 ij 

and k.,k. lie in the same sextuple. Similarly, when i = j , |k. + k.| = 2n_ 

I X J u 

so that the have a common value, C, i - 1,2,...,3N. Thus 


(6.4) 


A. . = 
ij 


C, if i = j 

A, 


if i 4 j and k, ,k. c T , n = 1,2,...,N, 
~i n 


Tt follows from (4.10) and (A. 37) that all A.. >0, hence A > 0; further- 

ij - - 

more, hypothesis (H c ) is equivalent to C > 0. 

It is also possible to determine other relationships among the A„ when 
l<_ ,k_. lie in different triples. Such relationships are not required to study 
the classical hexagonal cells but are given in (6.6b) below when N = 2. 


Remark 6.2 . In the context of (N,Mq)= (2,28) in Remark 6.1 there are nine 
distinct positive values of |k.+k.| for i,j e {±1 ,±2 , . . . ,±6} . Therefore, there are at 
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n. 


most nim distinct positive a^. “ a ^ . . One f inds that * n^ * C, 

0 1 

i = I , 2 , . . . , A , and 



'Mi 

“ a, 34 - 

a 56 ’ 

a l-6 

85 a 23 

" a 45 ’ a 14 * 

a 36 * 

n 2-5 * 

(n.6) 

<! 1 - 4 

a 25 * 

a 3-6 

» a 16 

" a 2-3 

= a 4-5 * a l-2 = 
n 2-6 

a 3-4 = ' 

a 5-6 ’ 

a 1 $ 

= a l-5 * 

a 35 

" a 24 = 

a 46 = 




‘l ) 

,l l- 3 * 

a 26 

* a 2-4 

* ft V5 

** iii /• ) ami 

‘4-6 

1-.1 

“ a i.i ‘ 

? f to! 

1 i i t hat ! 

he A n 

satisfy 






fa) A j j 

A 1 5 

a ;*4 

‘ A 2fi “ 

A '35 = 

\c, a 




o») a 12 

“ n 34 " 

A 56 * 

A 14 - 

A 25 = 

A 36 ' A 16 

A 23 = A 

45 * 

The fi 

•lat 1 onshipt 

> (6.6b) 

are 

needed 

for a 

complete analysis of 

"exotic 

So lilt 

ions when i 

i » 2. 







1 

I'rom (5 . 4h ) 

and (6. 

4) we 

get c 

■0 i) * 

c(iO + d(P), 

where 



m, . 7 ) 


3N 


l 5 *,! 


l( ' 1 K !. C " ** + \ A 1 ^i^i+N + " A+2N + l ‘ , i+N P 'i4-2N ) 

i-l 


, 2 n 2 


') ? 


and cl ( ) denotes the contribution to the sum in (5.4b) of terms . rfi-' , 

for which kj and lie in different triples. Note that d(b) > 0, 

. K 3N , since A.. >0. Thus, £(3) > f(S), where f (S) is defined 
ij " ~ 


in (5,6) and 


( 6 . 8 ) 


f(3) 


(q(6) + c<6))/|3| , if 3^0 
0 , if MO 


The functional f and its critical points play a key role in the determina- 
tion of stable, subcritical hexagonal solutions. Since (H^) is equivalent 
to C > 0 in (6.4), the functional c also satisfies (H c ) , so that 
i.emma 5.1 is applicable to both f and f. 
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Lemma 6.1, The nontrivial critical points, 3, of f satisfy 
2 2 2 

! n * ^‘n+N “ * n+2N * n * 1>2,..,,N. Moreover, f assumes its absolute 
minimum, 


(6.9) 


f Q “B“/9C l with C } “ C + 2A , 


2 o 0 

at those critical points for which all nonzero 3 satisfv P » 4B /9G7. 

n * n i 


P r o o f > At a nontrivial critical point we have (see (5.8), (5.9)) 


( 6 . 10 ) 


0 - \('\ 2 ; jjf;(P) - -2f(3)H i + x « 1,2,..., IN. 


let T' n * (k n ,l< n+N ,k n+0N ) be any triple and let i,.j,m be the indices 
(n,n + N',n + 2N) written in any order. Multiply the i— equation in 


(6.3 0) by L' ho j 


th 


equation by and subtract to get 


( 6 . 11 ) 


0 - (i 2 ± - (ij)[-2f(P) 4* C(fjJ + 0*) + A p.2 ] , 


By making use of the equivalence of (i) and (ii) of Lemma 5.1 applied to 


one sees as in (5.14) that f(3) < 0. Hence (6.11) and (H^) imply 


L hat 
have (• 


2 

1 

>2 

n 


4 

2 

3 ... - 0‘\ OM , n = 1,2,...,N. Observe that if k. <: T the 
n+N n+2N’ * ’ * « x n 


Since n and the order of i,.j,m are arbitrary, we 

,2 


.th 


equation in (6.10) involves only $ n , ^ n +N anc * *'n+2N • Since we 


may change the signs of any pair of these three 's without changing 

tii ** 

tlie i — equation, we may suppose at a critical point of f that 


t'n + 9 N ~ 3 n+N - $ n > n = 1,2,...,N. Then the three equations in (6.10) 

corresponding to each T (1 become identical and (6.10) reduces to N 

equations for 3^, n - 1,2,...,N. We suppose that exactly Mq of the 

3 are nonzero and reorder the indices so that 0 ^ 0 if n - 1,2,..., ML, and 
n n 0 


3 =0 if n = M + 1, 

n 0 


,N. Then (6.10) may be replaced by 
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((>.12) 0 - -2f(B) + B0 ± + Cj^, 

whore C.j » C + 2A. When tfy" 1 one solves (6.8) and (6.12) to obtain 
t >j » 4B 2 /9 Cj, f(3) * f Q r ~B^/9Cj, . When Mq ^ 2 one subtracts the equation 

tor ;.j from that for to gel the. MqO^- 1) /2 equations 


(6.13) 0 - (i^- i3 ) [B + C 1 (B ± + ftj)l » J - i + 1,... > M 0 ;1 * 1 » • • • , M Q * 

It is easy to deduce from (6.13) that the Bj's either are all equal or 

assume exactly two distinct values. When the B.'s are all equal, the system 

J 

((>.8), (6.12) becomes a pair of equations for 3^, f(B) and one finds 
) o ' } 

that f-J “ 43 /9C“, f(fci) ■ f Q . In the case of exactly two distinct ;■ , 

we suppose hj ^ 3*2 with p^ of the Bj's equal to ,* ; j and p 2 of 

the i-'.'s equal to 3 2 > \\ + P 2 = m q. Then the s y stem (6.8), (6.12) reduces to 


(6.14) (a) f(3) = £B( Pl i^ + p 2i i 2 ) + fc i (p l aj + P 2 4^/ 3(p l B l + p 2h ) 


(b) 2f(;0 - + C 1 p 1 , i — 1,2, 


Since , B 2 are different and nonzero we seek a solution iri the form 
= S3,, 8 f 0,1. Using (6.14) to express f(3) and p in terms of 
s, one finds that = -b/c.j (1 + s) , f(3) - + s) 2 and the 

solutions are determined by the roots, s, of 

4 3 

0 = -p 2 s + 2p 2 s + 2p 1 s - . 

The latter equation has exactly two real roots s^, s 2> which satisfy 
U _ Si |, 2 < s 2> If is the value of f corresponding to 

s^, i = 1,2, then one shows that f ± > f Q so that these solutions do 
not give the absolute minimum of f. 

3N 

Recall that f(3) f(8) for a H 8 efi and, in addition, observe 

that f(B) = f(6) = f Q = -B 2 /9C 1 at points, 3» of the form 
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( 0 . 13 ) 


Since 
of i’, 


» 0 if i 4 N, n + N, n + 2N 
, ’ n " t: n+N “ l ii+2N " 2B/3C j » 11 " * 

is the absolute minimum of f, it is also the absolute minimum 

i.u. 


(0,10) i’ n - min f O’) . 

u i#. ir lN 

it follows from (6.16) that each point, of the form (6.15) is a critical 

point of f. Moreover, since Cj > 0, one can show that the matrix E in (5.13) 

6 N 

is nonsingular at these points; in fact, det E (fi^C^/2) . Thus, according 
to Theorem 5.3 each of the points (6.15) generates a subcritic.nl solution, 
v = v(n,N), of (*),.. stable in S , Note that because of (6,16) and Lemma 3.1, 
there are no other solutions in 5^ generated by solutions (f:,") of (5., 5) 
with *" < Tq 5 2f^, 

Re mark 6.3, One can, of course, also consider the solutions v(n,N) 
as solutions of (*) in H. The stability of the v(n,N) in H is de- 
termined to lowest order by the eigenvalues of the 6N * 6N Jacobian 
maLrix of the full selection equations (3.24) at y - r - 0. One 
finds ns in [T, pp. 642-643] that all but two of these critical eigenvalues 
are positive and, because of the invariance of the equations (2.1) under 
translations of Lhe (x,y)~plane, the remaining two are 0. Thus, the 
stability arguments in [1] apply also to the hexagonal solution v(n,N). 


We shall call a solution, v, of (*) a hexagonal cellular solution 
if the leading term in v has zero component across the vertical faces of 
a right hexagonal cylinder Z and also across the vertical faces of cells 
obtained from Z by repeated reflection across the vertical faces ( the 


,* * ♦ 
- * 
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UK i h of 2 is parallel to the x-axis and the crons sediomi ?. * ? () » 

t j 

~ • > y, n ,1 ') are regular hexagons). For example, the solution v(n,N) 
generated by (6.11) Is a hexagonal solution (note the shape of the 
streamlines in (10, Fig. 1]; see also [3, 116]). One can show (e.g., 
see | j, H6|) Lli.it ij* *■ 1- tp + tp has xero component, across the 

vertical laics of i whose crons section x - 0 is the hexagon with 
center at (x,y) * (0,0) and vertices at g(4r/3 ^)k^, * T n - Clearly, 

tlte same Is Lrue of tp, corresponding to k^ < *T , hence of ’ * T + v ■ 
Furthermore, the. flow r , has the positive x-direetion along the x-axis. 
Thus, we see that the leading term in v(n,N) has this hexagonal structure 
and, since 0, '■ 0, th e flow I s upward along the 2 - axis when 

, 0 and downwa_n(_ when r ■' 0. 

one may also investigate the existence of exotic solutions in S t 
for general N by the methods of the present section. To determine the 
suability of exotic solutions, however, requires Lhu verif ication of cer- 
tain inequalities among, the coefficients of the functional f in (5.6ul. 
This is Illustrated in the following remark for the case N = 2. 


Remark 6,4 . Besides the simple hexagonal solutions determined above, 
one obtains in the case N * 2 additional solutions corresponding to 

(6.17) (a) * * -2B 2 /9(C ] + A^ = -2B/3(C X + , i - 1 6 

(b) i = + (i + + s 2 )Aj](5 x ,f3i = ~ -B/Cl+s^) (C^-A^) , 

= ^4 = “ s i^i» 

whore - A^ 2 4- A^ + A X( . (see (6.6)). Here s^, 0 < < 1, is a root of 

(6.18) 0 » 2(C X - A l )(s 3 + s 2 + s) - (C + 2Aj) (s 2 + l) 2 . 


One finds that Lhe existence of s^, lienee of (6.1.7b), as well as the 
stability of both solutions in (6.17) depends on Lhe sign of - 7A^ , 
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I'lu* solution (6,17a) corresponds to Liu* l irst exotic solution tease da) 
in [ 10J . 

Observe that In unch of the solutions (6.15)* all h'^'a corres- 

ponding to a given triple, T n » are equal, The functional f, however, 
lines not change If we change the signs of any two o^'s corresponding to 
the same triple. Thus, each of the hexagonal solutions generated 

by (6.15) yields three additional hexagonal solutions. One can show that 
tin* lour solutions obtained this way are translations of one another. 

Moreover, aU of the solutions, v(n,N), generated by (6.15) (for 
u *• l, N; N in a suitable, unbounded sequence) are, at least to 
first order, rotations of v(l,l). 

our main results for classical hexagonal cellular solutions are sum- 
marixud In t he next paragraph and hold under the hypotheses C 0, B i 0 
(see Remark 5.1). These hypotheses are independent of N and are analogous 
to Lhu minimum hypotheses required for a bifurcation analysis at '• when N * 1. 

Hexagona l cellular solutions . For each N in a suitable unbounded 
sequence there are 4N solutions of (*) generated by absolute minima of 
the selection functional, f, These solutions are subcritical and stable 
In S (| * S. n (N). Each of these solutions exhibits the classical hexagonal 
cellular form with sice independent of N. The stability of, e.g., v(l,N) 
in S, (N) shows that the hexagonal cellular solutions are, in particular, 
stable to perturbations in "directions” corresponding to N critical wave 
vectors. Thus, letting N range over the unbounded sequence, we obtain, 
in a sense, the stability of the classical hexagonal cells in infinitely 
many such critical directions. 



There 1 b no attempt. In the present paper 


/ . (Ii nic ludijijt jenurr R k . 
to obtain "nil oi the local so lu Lions near \ ■* i.j" of the Benin <1 
problem with symmetric boundary conditions even in the simplest of i,n.e».. 
Jhe motivation has been rather to provide a first step toward showing that 
Lhe hexagonal cellular solutions are the "preferred" subcrltie.nl solutions 
oi the Benard problem in physical situations with temperature dependent 
material properties. In fact, the recent results of Buzano and fiolubit.sk-.* 
(2) and Uoiuhitsky, Swift and Knobloch [ 5] indicate how difficult it would 
bo to obtain "all of the local solutions near ' * even in the case 

in Section (i when dim M * 1,1* fn [2], [5] those authors consider situa- 
tions corresponding here to the case in Section (> of one triple of critical 
wave vectors, i.u., dim M * 6 and, by an application of group theory 
and, in [2], also singularity theory, they obtain "all of the local solu- 
tions'' of a six-dimensional problem P. (One assumes that V corresponds 
to the f IniLe-dimunsionai problem generated from the Benard problem bv 
means of the Lyapunov-Schmidt method relative to the first eigenvalue of 
the linearized problem.) The detailed results in i 2 ] are of particular 
Interest because they show for the Benard problem that the mathematical 
possibility exists of having stable subcritlcal hexagonal- type solutions, 
stable supercritical roll-type solutions, and a Lhlrd type of solution 
that provides a transition between rolls and hexagons. There are, of 
course, some difficulties encountered in carrying over the finite-dimen- 
sional results in [2], [b] to an infinite-dimensional mathematical model 
and many such difficulties and their interpretations for the Benard 
problem are discussed in [2, 511], The most pertinent such difficulty 
relative to the method presented here is the fact that the detailed nature 
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of tin* rt.*Hu] ts In [2], [5] are highly dependent upon thy relatively low 
dimension of the problem V whereas the basic results of Busse [1] are 
essentially independent of the dimension of any underlying f inlte-dimen- 
sionnl problem. One of the main goals of our study of the Benard problem was 

to develop a rigorous stability method useful In a setting that also is 
independent of the dimension of any underlying finite-dimensional problem, 
flu* results of Section 6 show that this goal has been achieved and that 
in our approach Lhe selection of stable suberitical hexagonal cellular 
••■■luiions i-. closely related to a minimization condition on the general irn-d 
dissipation. As in earlier work on the Benard problem (e.g,, see [1; 

I)];, there remains in the case of temperature-dependent material proper- 
t ies Lhe problems of finding a strict physical interpretation of the 
generalized dissipation and a description of the actual selection mechanism. 
Finally, we note that the methods introduced here can be modified to 
vield also the description of stable supercritical states and the stability 
1 1 * I at hmships between roll-type solutions and hexagonal cellular solutions, 
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AEESEL^* Here we justify equation (2.23) and prove Lemma 3.1. 

Fir«L we allow that tin* eigenfunctions {'; nq ^} in (2.20) can he sealed 
with constants independent of ,j so that (2,23) holds. In fact, we may 
assume that each <J pq lias been scaled by a constant depending only on p 
and q such that 


(A.i) 


2 

~2 

i 

P 


'U 2 
• - 1/2 


[D 2 ^ 


]“ dz 


fX l U 2 






d 2 

where I) “ - v - ) . (The integrand on tlie left, in (A.I) is not zero, bv 

dz p 

uniqueness of the initial-value problem 1)^ = 0, if>(i) - >|i'(|) = () (see (2.12e)),) 
From (2,3) and (2.21) we get 


(A. >) 


(d' Pqi ,'!») = j(p,q,.j ;r,s,t) 


< n/, h f 2 ’ ,/,J 2 i( k . - k ) . 


L 


0 - pj *vdydx 


4' r T ** 

’ ‘ S pr 6 j C J Cp,q,j ; r 

Her*-, since A pC ^ = d‘ Pq and 4*/^ - <|> Pq are real and independent of j, 


(A. 5) 


J Cp ,q , J jr,s,t) 


4 4 

1/2 ' 

'- 1/2 

3 
Z 

m=l 


(k . • k ) a ^ rstr 
- PJ - rt m=1 m m 

1 

+ T -i-^pqj Jl I rst 4-1-rOc k )* p V 8 + ^ pt| V s Ud 
dz*\n dz + Pr Kk P j k rt' 0 4 ’4 + dz T 4 dz'4 J | d 


From (A. 2) we see that J is needed only when r - p and = j , Then we may 
integrate by parts in (A. 3) making use of (2.12) to show that 

rl/2 

d(p,q,j;p,s,j) = ( i pq j (^ q <J> 3 S + <^ q '-i>J S )dz . 


(A. 4) 


Since both ij) pq and (fi ps satisfy (2.12) we have, after integrating by parts, 


0 = ()J - 1J )o 


rl/2 


pS pq P ■'-1/2 


(c|) pq (f) pfa + <J) Pq (f) PS )dz . 


(A. r >) 
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*pqj 


I’huK (A. 3) shown that if u ? ii (l.e., if s # q) then y' u and 

pa pq 


, PkJ 


are orthogonal in the sense that 


(A.O) 

Jn particular, 
(A. 7) 


1 1/‘ 


U’ 


op*?* + xi* • 


h 4 3 


3 4 


•J(p,q,j;p,s,j) * ft J (p.q,J;p.q,j). 


hut irom (2.12a), integration by p<krts and (A.l) we get 


(A, 8) 


A ( p » q » j ; p » q * j ) ■ 2o~ 2 f . 


p ’'—1/2 w " 4 ti‘ 

Combining (A. 2), (A. 7) and (A. 8), we obtain (2.23). 

Next we give a proof of Lemma 3.1; some aspects of the work is closely 
related to corresponding steps in [1] or in [10]. According to Lemma 2.1 
the operators l. and M are bounded on H. Since L is also compact, it 
is easy to see that K is bounded on M X . If v • H has the form (3.3) 
and A is any bounded linear operator on H, then Av may be computed 
Lerm by term in the sum so that the formulas (3.7) follow easily. The 
positivity and self-adjointness of K are simple consequences of (2.23), 
(3.7) and the fact that the p are real. 

pq 

Part (ii) of Lemma 3.1 follows easily from the definition of M if 
we show that 


(A. 9) 


.y 


0 = (M/) m ,ip) , for |m| , 1 j | - . 

Po 1 


Since ip^ and <|> ^ ~ are even functions of z (see [10]), (A. 9) 

follows from 


(!#'\V J ) = 2 


— 6(k + k.) 
^ r 4 r 3 a 1 « 2 -m 


1/2 


z<J>. <j> 0 dz = 0 . 


- rj - i /2 '«' 3 

The assertions (iii) of the lemma are obtained from (3.2), (3.3) and 


the identity 
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43 , 


J (u*Vv)*w 


J {V*fu(v*w)] 


(u*Vw) *v) 



(u’Vw) • v 


This Inst identity is easily verified for smooth u,v,w <• H and is proved 
in general by a standard limiting argument using the boundedness, in u,v,w, 
of tlie functional (!(u,v),w): 


(A. 10) 


OHu.v) * w ) | 5 const || u U || v || |jw||. 


(The inequality (A. 10) follows from (3.2) by application of the Schwarz and 
Poincare inequalities.) In addition, one may use (A. 10) to show that all 
assertions in part (iv) of Lemma 3.1 are consequences of (A. 11) and (A. 13), 
oeiow. 

On substituting (2.20) and (2.21) into (3.2) we are led to (recall that 


'0 


and p„l are suppressed: a - a, k . * k,, etc.) 


P0 


‘P o J 


(A. 11) Ct'(!> Pqm ,tp j ),^ n ) - - 


4tt 

u l a 2 


i$(k pm + k j + k n )I(p,q,m;p 0 ,j,n) 


where 


(A. 12) 


i(p,q»m;p 0 *j,n) 


rl/2 


-9 d< ^ q ■? 

' - 0 ~ l (k *k ) — 7- — (<r*$ u ) 

^-1/2 ( P ~P m dz 


+ <J> 


pq 


• , d <J>~ d(|) d(J)„ 

a 4 (k-k)— J J 1 J 
-J ~n dz 


2 dz + dz ^3 + 


d ^4 "1 . 

dF +*| j dz 


The right hand side of (A. 11) is zero because of the 6 term, except when 

k + k. + k = 0. In this exceptional case the vectors k., k and k 
pm j ~n -3 ~n -pm 

I orm an isosceles triangle so that k *k. * k *k . Consequently m and 

-pm -j -pm ~n 

j may be interchanged in (A. 12) without changing I. In this case then, 
from (A. 11), and part (iii) of Lemma 3.1 we have 
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44 . 


To prove the formulas of part <v) of the lemma we take •/ * 
and calculate the various terms. Now 

(A. 14) » 2 f z^ pqm dz - 5 6, b n . 

Jo 43 PPn 1"' 0< 1 


h l-i 


ilere 

„ 2 f 1/2 p n l p_q 

IA - ,r ° b o, Li/i^ *3 dz ’ lq l “ 1 ’ 2 ”-- 

is real and = k()(_^) = 0 since and are even. Then 


N 


Mtf» - F. fi. 

IJI-I' 1 


Z b n i/f 

|q|-» 0<I 


(A. 16) 

Similarly, from (2.24) we are led to 

_ n 00 * -Pn^ 0 

(A. 17) H*il> n = 2 b* U 

I q 1=2 0 « 


Pn*1.1 


where 


(A. 18) 


* _ 8 ir 


,1/2 P„q p p l 


°q ~ <\n 2 J 


• 1/2 


ztp^ cj> 3 dz, | q | = 1,2,... 


is real. Since M : r c ill , KMijj may be obtained from (A. 16) and (3.7): 


(A. 19) 


whe re b 


N 


KMiji = l 6. 2 b n tj) 

i j 1=1 J |q|-« ° q 


Pn9H 


Oq 


V_ _(y_ _ - H-.) ^b- . From (A. 17) and (A. 19) we have 
Pnl PnT 1 0f l 


0 p 0 


(A. 20) 


-n N f 00 ~ * P 0 m i _P 0 qn 1 

(MKMi|j,^ n ) - E 3, 2 b n b’ (tp U «bJ3 

Li l = i tl m l » I q 1 ~2 0m0q J 0 - n 


Here 


(a. 21) b = E b b“ 

0 jq J=2 0q 0q 


is real . 


This proves (3.8). 


M Z 


Wo shall require 
(A. 22) 

where, hv (2.21), 

(A. 


omoiim. i 
of poor c. 


( , K*; 1 , i m ),l pc]n ) * £(k. + k n + k pn )T 2 (P 0 ,i > m »P*q»n) , 


45 . 


, n r / • -v - 4ir 2 T -2,. , N d,! 3,ro „pqn , fdf m ,nqn) \ , 

23) i » in * p * q » n ) ^rVar 1 • “^,dr* r • 

Wi* are interested in wily when the vectors k., k m> k form an isosceles 

triangle, otherwise the e-factor in (A. 22) is zero. From this trianrle we 

12 12 2 

see that k.*k = k *k = - v a and k,*k = •„ i - . In this case, 

j pn m pn 2 p j m 2 p 

(A. 23) leads to 1 2 (p Q » j »"),p PI >n) = ^Pq’P*^ > where 


dja 


(A. 24) l 3 (p 0 ,p,q) - - 


4tt 2 , 

(1/2 f 
y | 

r 2 'i 

•i o 

1 _ d __E. 

th 

' 4 l' x 2 J 

-1/2 

2 ^2 

k ‘ > 

dz 


x^.lxV * 


pq 


2o 


2 dz dz '3*3 


4- ■*, • i 


.pq 


+ -:p 


, d 2 fj).> d>j. Pq d(K d}., ] 

1L + J.,1 (] ,pq + JA ( hPq 

, ,2 .2 dz dz q 3 + dz P 4 

2o dz ! 

- 1 


dz 


is real and depends on j and m only through = jk^ + k^J. On the 
basis of (A. 22) and (A. 24) we have 

(.A. 23) ,i|™) ,i[' Pqn ) * iS(k + k + k )I~(p fi ,p,q) . 

~ j -m -pn j u 

From (3.6), (3.7) and (A. 17) we have 

00 -jl p n qn 

(A. 26) (MKF(i»,l") = (F(«),KM*)") - Z b* CFO,) ,0 ), 

!q 1=2 uq 

where b" r u (p - y..) is real. Furthermore, from (A. 25) with 

Oq P Q q P Q q H 1 Oq 

p = p 0 we get 


- p o qn 


) - S B 3 m 6(k + k m + k n )I 3 (P 0 .P 0 .q) • 
| j | , Jm| = l 


(A. 27) 
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Combining (A. 26) with (A. 27) we obtain (3.9) with real hj given by 


(A.2M) 

b 1 " | j p . t} b 0q I 3 (p 0‘ P 0 ,fI ^ * 


In a similar 

manner we may utilize (A. 19), (A. 25) 

and part (ili) 

1. omnia 3.1 to 

establish (3.10) with real constant 

b 0 given by 

(A. 2°‘ 

IX) 

b 2 * | q |_ 2 1> 0 q I 3 <l ’0 > ' , 0’ < l) • 


l.tjiiation 

(3.11) follows easily from the observation that 

(A. ID) 

( p ( ; , i ,m ) fi ,n ) ■= ( k . + k t)) + k (i )b j , 


where 



(A. 31.) 

^ rf 2 fl/2 2 

>’3 - ^ [ 1/2 ( 'V V > ° • 



Next we consider part (vl) of the lemma. From (Hi) of Lemma 3.1 
and the bilinearity of <!> we have 


(A. 32) 


,KF(i|j) ) ,{jj n ) = ~ (^( 4 , ,KF< F) ) 


We may obtain the last inner product by means of Parseval's equation as 

follows. From (A. 25) the coefficient of in the Fourier expansion 

of ‘Hr' 1 , !’ 11 ) is A(k. + k f k . )I„(p_,p,q) , while the coefficient of 

i ~ n *v p n j u 


pqh . 


(A. 33) 


in the Fourier expansion of ,t|j m ) is 

,S(k i + “!Spl, >|, pq <|J p<I “ l , l > " 1 L <P 0 ' |, O)- if <P.1> * V> 

0, if (p,q) = (p 0 ,l) . 


Consequently, 
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_ l( , _ t | M . | p X ” ( j? p ) 

(A ' W - E 0 «a j+ k n+ k ph )6(k 1+ k m -k pli ) 

where Eq denotes summation over (p,q,h) in the same set of integer triples 
as in (3.7). 


Oiven j and n, the only way a term in the sum on the right in (A. 34) 
can be nonzero is for 

(A. 35) k , * -k. - k * k. + k . 

pli -j **» x *viii 

These relations determine p and h completely, in terms of p^.j and n 

(or In terms of Pq, 1 and "0 so that only q need be summed in (A. 34). 

The relations (A. 35) also require that either k. = -k. and k = -k (i.e., 

«< x ■•< j - m * n 

i «■ - j and m = -n) or k . * -k and k = -k. (I.e., 1 = -n and 

"X *• ti -m 

n. * =j). It follows that 


(A. 36) 6(k.+k + k ,)<$(k.+k -k ,) = 6(i + j)6(n + m) + 6(i + n) 1 (m + *•) 

~] -n ~ph -an ~ph J 

- 6(1+ j)5(i + n)6(i ~ m) 

If we combine (A. 32), (A. 34) and (A. 36), then we obtain (3.12) with non- 
negative constants a , given by 

P Q jn 


00 


(A. 37) 

V" 

1 ,J pq* ,J pq“ 1 V I 3 (p 0 ,P 
q=q x 

*q) * 

i .1 1 » 1 n | = 1,2,... 

where q^ = 

1, if 

P i 1 P Q and q x = 2 if p 

ii 

-o 

o 

Note that a 

P 0 3 n 

depends on 

j and 

n only through p, i.e., 

through a - |k, + k 
P '■! ' n 


In particular we have 

( A • i 8 ) a « = a , — a * • \ * \ • 

P 0 in p 0 n.i p Q (-j)(-n) 


Furthermore , when .j = -n we may, for convenience, define 
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4H. 


(A, ')')) 


a P 0 iH) “ 


(Thu sum in (A. 37) is meaningless in tills case, since (2.12) lias no turn- 
tilvlal solutions when n * 0 (* jkj I- k^_^j) so that 0 4 for 
any p) . Because ( p » q ) ^(p^,!) in (A. 37), we see when n 4 -j that 
a^ * 0 if and only if 1 g (P q ,P i f l) “ 0 for all integers q with q ' 


A i ' 1 <j u ) w 1 ini g e me at . The authors wish to thank Professor E, A. Connors for 
acquaint inn them with the number-theoretic results used at the beginning 
of Beet Ion 6 . 
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